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We study the cosmology of a toy modified theory of gravity in which gravity 
shuts off at short distances, as in the fat graviton scenario of Sundrum. In the 
weak-field limit, the theory is perturbatively local, ghost-free and unitary, although 
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CN ■ likely suffers from non-perturbative instabilities. We derive novel self-inflationary 

<-> ■ 

solutions from the vacuum equations of the theory, without invoking scalar fields or 
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other forms of stress energy. The modified perturbation equation expressed in terms 
of the Newtonian potential closely resembles its counterpart for inflaton fluctuations. 
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. The resulting scalar spectrum is therefore slightly red, akin to the simplest scalar- 

m ', 

driven inflationary models. A key difference, however, is that the gravitational wave 
spectrum is generically not scale invariant. In particular the tensor spectrum can 
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have a blue tilt, a distinguishing feature from standard inflation. 
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What if gravity becomes weaker — and perhaps even shuts off — at short distances? 
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This tantalizing possibility should have dramatic implications for early-universe cosmology. 
An immediate question is whether high energy inflation can still take place — after all, the 
accelerated expansion of inflationary cosmology relies on the backreaction of vacuum energy 
on the geometry. But maybe shutting off gravity at high energy can obviate the need for 
inflation and offer interesting alternatives. If no gravity means no cosmic expansion, the 
universe could conceivably start out in a quiescent state, having plenty of time to thermalize 
and homogenize. 

The most widely studied classical modifications of Einstein gravity, on the other hand, 
such as those arising in theories with extra dimensions, all have in common that gravity is 
stronger at short distances. This can be understood in a variety of ways. First, stronger 
gravity is seen in the Newtonian potential becoming steeper than 1/r above the compacti- 
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fication scale: 
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Stronger gravity can equivalently be understood at the level of the 4d effective theory: 
the various moduli describing the size and shape of extra dimensions mediate attractive 
yukawa forces, thereby enhancing gravity. Finally, stronger gravity leaves an imprint on 
the cosmological evolution. For instance the modified Friedmann equation [1, 2] in the 
Randall-Sundrum brane- world scenario [3], 
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implies faster expansion than in standard cosmology. But this is precisely what we ex- 
pect from stronger gravity: for a spatially-flat universe to expand forever, as it must, the 
expansion rate must be higher to overcome the stronger gravitational pull. 

This paper considers instead the implications of weaker or fading gravity on early-universe 
dynamics. The idea of "asymptotic safety" for gravity [4] has been exploited in various 
contexts [5]. It has been proposed to resolve black hole singularities, either through non- 
perturbative a 1 effects in string theory [6], or in models inspired by the non-local form 
of the open string field theory action [7-9]. Weaker gravity is featured in the recent "fat 
graviton" [10] proposal to address the cosmological constant problem. An interesting search 
for fading gravity in brane- world constructions is described in [11]. 

We focus on a particular toy model of fading gravity, where the graviton propagator shuts 
off analytically in the UV: 

Here Q goes to zero on scales much smaller than L and becomes trivial on large scales. 
(Alternatively one can suppress gravity with the addition of a massive or Proca spin-1 field 
that couples to matter with charge m/M P1 [12]. As in electromagnetism, a vector- mediated 
force is repulsive for like charges, quenching the effective gravity. See [13] for cosmological 
applications.) In order to avoid adding new poles in the propagator, which would invariably 
introduce ghosts, Q must be analytic, e.g., 

g (DL 2 ) = exp (— D 2 L 4 ) . (4) 

More generally, it can take the form exp (/ (QL 2 )), where / is an arbitrary analytic function 
consistent with the asymptotic behavior of the form factor. This mechanism for shutting off 
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gravity at short distances is inspired by the "fat graviton" proposal [10]. Analytic propaga- 
tors of this form appear in Polchinski's formulation of the renormalization-group flow [14], 
and have been considered to regularize divergences in gauge theories [15]. In the cosmolog- 
ical context, modified propagators have been proposed to derive bouncing cosmologies [16]. 
With normalized kinetic term, the action is reminiscent of the tachyon sector of bosonic 
open string field theory as well as the scalar field action in p-adic string theory [17]. 

The spatial and temporal nonlocality of our form factor leads to an action with infinitely 
many time derivatives, which therefore lacks a well-defined initial value formulation and 
likely suffers from non-perturbative Ostrogradskian instability [18, 19]. Nevertheless the 
theory seems remarkably well-behaved in most regimes of interest for this work. 

In the weak field limit, for instance, the theory is perturbatively ghost-free, unitary 
and has a well defined initial value formulation. Indeed, despite arising from equations of 
motion with infinitely many time derivatives, all perturbative solutions are specified by a 
finite amount of initial data. This can be understood by the localization procedure [19]: 
any analytic non-local theory with local unperturbed lagrangian can be recast through field 
redefinitions as an action that is local in time and reproduces all perturbative solutions. 
Since any curved background is locally flat, in this sense the UV limit of our theory seems 
under control. 

Even perturbatively, however, the theory is not without worries. It is well-known that 
analytic modifications to the propagator lead to acausality. At the classical level this is seen 
directly in perturbative solutions having support outside the past light-cone. In quantum 
field theory, the culprit for acausality is the essential singularity of (4) at infinity in the 
complex frequency plane which spoils the analytic properties of scattering amplitudes [20]. 
Be that as it may, at a philosophical level acausality at high energy is arguably a mild sin. 
It may even be a blessing for early-universe dynamics, in particular for the homogeneity and 
horizon problems of standard big bang cosmology. 

A. Self-Inflation 

The modified Einstein equations, whose weak field limit describes a graviton with prop- 
agator (3), are of the form 

g- 1 (DL 2 ) G a(3 + O (R 2 ) = 8nG N T af3 , (5) 
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where the curvature-squared terms ensure the Bianchi identity. Focusing on the leading 
term, the form factor Q is recognized as an effective Newton's constant, 



which indeed vanishes at short distances. 

The corresponding modified Friedmann equation offers us a pleasant surprise: an accel- 
erating universe without invoking scalar fields or any other form of stress energy. This novel 
inflationary solution instead arises from the vacuum equations of the modified gravity ac- 
tion — a self-inflating solution. This is interesting phenomenologically since a lot is known 
about inflation with scalar fields: the required fine-tuning on scalar potentials, their generic 
predictions [21, 22], etc. Our theory, however, is a genuine spin-2 modification and there- 
fore cannot be rewritten as Einstein gravity plus scalars. Thus the resulting accelerating 
solution is different than scalar-driven inflation and, as we will see, leads to distinguishing 
observational signatures. 

To understand how self-inflation is possible, a key observation is that de Sitter space, or 
any Einstein space for that matter, identically satisfies OG'Q = 0. Thus the action of Q^ 1 
on the Einstein tensor is trivial, Q~ x (OL 2 ) G^ = G^, no matter how small the dS radius 
is compared to L. 

Now consider a geometry which is nearly de Sitter, but slowly-evolving, in the sense that 
\H/H 2 \ <C 1. That is, our ansatz is R a p = 3H 2 g a /3 + r at3) with r a small correction of 
order H. Since r a p is slowly-changing, we can neglect its time- derivatives, which amounts 
to dropping terms of order H, H 2 , etc. Schematically, we find 



Thus, although r aj3 constitutes a small deviation from exact de Sitter, its contribution to 
the equation of motion is "boosted" by a factor of Q^ 1 (8H 2 L 2 ) 3> 1 for HL > 1. 

Our solution therefore can be understood as arising from a compensation between the de 
Sitter H 2 term and the "boosted" H correction: 




(□) = g (DL 2 ) G 



N , 



(6) 



Q- 1 {UL 2 ) r a , ~ Q- 1 (8H 2 L 2 ) H . 



(7) 



3H 2 + Kg- 1 (8H 2 L 2 ) 



(8) 



that is, 



H 



-G (8H 2 L 2 ) . 



(9) 
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The factor of Q (8H 2 L 2 ) on the right-hand side is indeed small provided that HL > 1, 
justifying the slow-evolution approximation. The latter can be checked explicitly by taking 
derivatives of (9) and noting that H 2 , H etc. are down by extra powers of Q. While this 
self-inflating solution relies on the higher derivative nature of the theory, nevertheless the 
approximate equation of motion (9) from which it derives is second order in the scale factor, 
just as in Einstein gravity. In this sense the dependence on the higher derivative structure 
is the most minimal possible. 

Intuitively, self-inflation can be understood as follows [44]. The fading of gravity arises 
from self-interactions of the graviton which become relevant on scales smaller than L. These 
self- interactions in turn give rise to the curvature-squared terms in (5), which manifestly act 
as a source for gravity. However, thanks to the filter function this source inevitably appears 
to gravity as approximately homogeneous in space and time on scales smaller than L — it 
therefore plays the role of an effective vacuum energy, causing the universe to inflate. 

Surprisingly, de Sitter space is an approximate solution of increasing accuracy for larger 
values of HL, and generic higher-curvature corrections to the effective action also become 
less and less relevant in this limit. The latter either involve derivatives acting on curvature 
tensors, which means they are subleading in H; or they are higher order in curvature, in 
which case they simply "renormalize" the form factor Q. This is in stark contrast with 
the usual situation in effective field theory where corrections to the effective action become 
increasingly relevant above the scale of new physics. This highly desirable property of our 
accelerating solution traces back to pure de Sitter being oblivious to the form factor. 

The inflationary era terminates when H ~ H 2 ; that is, when H ~ L^ 1 . The details of the 
transition to decelerated expansion have yet to be worked out through numerical analysis, 
which we leave for future work. Nevertheless, we know by construction that the cosmological 
equation reverts to the usual Friedmann law for HL <C 1. Moreover, when H ~ H 2 
the universe reheats to a temperature of order T ~ L^ 1 through gravitational particle 
production [23, 24]. Since there is no inflaton field that couples directly to matter, this is 
the only way to reheat. Therefore, we conjecture that the universe exits the inflationary 
phase smoothly and enters a radiation-dominated epoch governed by the usual Friedmann 
equation. 
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B. Perturbation Spectrum 

Our accelerating solution can be understood in a different way, which is useful for studying 
the generation of density perturbations. Anticipating an approximate de Sitter solution, 
consider the metric ansatz 

ds * = _ e 2Ht) dt 2 + e -2<Kt) e 2Ht dx 2 ^ (1Q) 

where H is constant and <fi is assumed small. At linear order, the latter is just the Newtonian 
potential. Our modified gravity equations can then be cast as an equation of motion for 0: 

3H<j) w 8H 2 Q (8H 2 L 2 ) . (11) 

In analogy with the slow-roll equation in standard inflation, 3H(p ps — V v , we see that 
the Newtonian potential plays the role of an inflaton, setting a clock for the cosmological 
evolution. (To make contact with (9), note that H = —H<j) ps —H<j>.) 

Density perturbations therefore originate from fluctuations in the Newtonian potential, 
&(x, t). In the limit where modes are well within the horizon, the geometry is approximately 
flat, and the equation for reduces to 

g- 1 (k 2 L 2 ) k 2 <$> k ^ . (12) 

Since Q is analytic, there are no new homogeneous solutions other than the usual plane 
waves in flat space, ~ e~ lkt / V2k. Meanwhile, in the opposite limit where modes are well 
outside the horizon, $ fc satisfies 

3H$ k ~ H 2 Q (8H 2 L 2 ) $ fe , (13) 

where we have dropped the bars for simplicity. Since Q (8H 2 L 2 ) 1, an approximate solu- 
tion is $fe ~ const. These two extreme limits coincide with the corresponding limiting equa- 
tions satisfied by inflaton perturbations in usual slow-roll inflation. By analogy, it follows 
that $ acquires a scale-invariant spectrum with amplitude A$fc ~ H/Mpi, corresponding to 
a curvature perturbation 

C~f, (14) 
where e = —H/H 2 is a generalization of the usual e parameter of slow-roll inflation. 
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Compared to the standard slow-roll result, Csiow-roii ~ H/e 1 ^ 2 , the amplitude of perturba- 
tions in our case is enhanced by a factor of e -1 / 2 , which traces back to the difference in the 
generation mechanism. In slow-roll inflation, perturbations in the inflaton field have ampli- 
tude {5(f) ~ H and are related to Newtonian potential through 5<p — — \j2/e + $ j . 
In our case, the Newtonian potential itself acquires (<&) ~ H, and this y/e mismatch trickles 
down to the final answer. 

The scalar power spectrum is generically slightly red in fading gravity. Indeed, e is 
proportional to Q(8H 2 L 2 ) and therefore increases with time, while H decreases with time. 
It follows from (14) that modes which exit the horizon earlier have a larger amplitude, 
corresponding to a red spectrum. More precisely, for the choice of filter (4), we find n s ~ 0.96 
on the observable range of scales, in good agreement with current data [25]. A blue scalar 
spectral tilt is also possible if one is willing to consider a non-monotonic Q. 

A key difference with standard inflation lies in is the spectrum of primordial gravitational 
waves: it is generically not scale invariant. This is not surprising since our fading gravity 
model is a spin-2 modification, and should therefore strongly affect the tensor modes. While 
it is nevertheless possible to achieve a scale-invariant tensor spectrum by judiciously choosing 
the form factor, in general this leads to a gravity wave amplitude on large scales well below 
the projected sensitivity of near future experiments. This follows from the fact that the 
tensor-to-scalar ratio, r, is proportional to e 2 in our model, as seen from (14), as opposed to 
e in standard inflation. Thus r is at best 1CT 4 — 1CT 3 in fading gravity. Hence, a significant 
level of gravity waves on large scales requires a relatively strong red tilt, which, if observed, 
would distinguish our model from scalar-driven inflation. The fiducial form factor (4), for 
instance, yields a strong red tilt. 

Another interesting feature is that the gravity wave spectrum is allowed to have instead 
a slight blue spectral tilt. Thus, if futuristic microwave background polarization or gravity 
wave detectors ever reached the 1CT 4 — 1CT 3 sensitivity level for r, observing such a small 
blue tilt would strongly hint at fading gravity inflation. This is because a slightly red tensor 
spectrum is a robust prediction of scalar-driven inflation: since H = — <p 2 /2 < 0, independent 
of the potential, the gravity wave amplitude ~ H/Mp\ is a decreasing function of time 
and hence a decreasing function of scale. Blue gravity waves also appear in the ekpyrotic 
scenario [26-28], but with such a strong spectral tilt that the large-scale amplitude is orders 
of magnitude below observable levels. A small blue tilt for tensors has also been obtained 
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recently from thermal fluctuations of strings near the Hagerdorn temperature [29], as well 
as in super-inflationary models [30] based on ghost condensation [31]. 

Perturbation analysis about the self-accelerating solution reveals a modified graviton 
propagator with 3 poles, one of which is ghost-like with a mass of order of the Hubble 
scale, indicating an instability of the solution. Intriguingly the ghost decouples in the limit 
HL — > oo in which de Sitter becomes an exact solution. To make the decay rate of the 
vacuum acceptably small, we must as usual invoke explicit Lorentz symmetry breaking at 
some high scale A LI [32]. (A Lorentz-invariant cut-off yields a divergent rate.) An explicit 
calculation for the form factor (4) gives A L i ^ 10 15 GeV. 

Once inflation ends and we enter the regime HL <C 1, the ghost further decouples and 
becomes a non-perturbative excitation. Indeed, by construction our theory has only two 
perturbative propagating degrees of freedom about flat space — the usual modes of the 
massless graviton. This decoupling is starkly different from theories of the form R + C\R 2 + 
c 2 RapR a(S [33] and F (R, R a pR a(3 , R aPl5 R a ^ & ) [34, 35] where the Weyl ghost persists all the 
way to the UV cutoff. 

The paper is organized as follows. We begin with a detailed exposition of the model 
in Sec. II. Given the complexity of the resulting equations of motion, we present three 
alternative derivations of the self-inflating solution. Exit from self-inflation and reheating 
are the subject of Sec. III. In Sec. IV we argue that our self-inflating solution is stable under 
generic corrections to the effective action. Section V deals with the generation of scalar 
density perturbations, their amplitude, and their spectral tilt. Tensor modes are discussed 
in Sec. VI. Section VII is devoted to a general discussion of pathologies of higher time 
derivative theories, in particular the Ostrogradski instability, acausality, and the breakdown 
of the initial value problem. We conclude in Sec. VIII with a brief discussion of future 
research avenues and open problems. 



II. FADING GRAVITY MODEL 

The cut-off graviton propagator (3) derives from the relativistic and covariant action 
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where □ = V a V a is the covariant d'Alambertian. Indeed, in the weak field limit, g a p = 
Vap + h a p, the Ricci tensor in Coulomb gauge is given by R a/3 = —Dh a/3 /2 + 0(h 2 ), and (15) 
reduces to 

^weak-field = ^ J d*X (h** ~ ^hj Q' 1 (DL 2 ) Dh a(3 , (16) 

where the propagator displays the correct massless spin-2 tensor structure but is cut-off in 
the ultraviolet. 

In order to avoid new poles, which would invariably include ghosts, the filter function Q 
must be of the form exp (/ (DL 2 )), where / is analytic. For concreteness, for most of the 
paper we will focus on the fiducial choice (4), Q (OL 2 ) = exp (— D 2 L 4 ), where the even power 
of □ ensures that filtering occurs independently of whether the given mode has time-like or 
space-like momentum. 

For this fiducial choice, the form factor in (15) can be expanded as (Q^ 1 — 1) /□ 
□L 4 + . . . and therefore does not include 1/D operators. It does, however, display an 
infinite number of time derivatives, which is a priori worrisome since temporally non-local 
theories generically suffer from instability and are bereft of an initial value problem. As we 
will argue in Sec. VII, however, things seem under control for all regimes of interest here, 
including the self-inflating solution discussed below. 

Putting non-locality aside, the action (15) nevertheless appears fine-tuned from the point 
of view of effective field theory. In general one expects O (R 3 ) terms, more general derivative 
structure, and so forth. We will come back to this issue in Sec. IV and show that our self- 
accelerating solution is in fact robust under the inclusion of a wide class of generic corrections 
to the effective action. 



A. Fading Einstein's Equations 

The derivation of the equation of motion from (15) is straightforward, modulo a few 
subtleties which are worth mentioning here. First, we must be careful in varying the form 
factor since □ and SO do not commute. For the fiducial choice (4) the final expression is 

G ' U 1 Sn(^—^-) -L 4 I' ds%^(SD-D + D-SD)g s -\ (17) 



□ j : : v : : j ./„ 
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where SD implicitly takes into account the rank of the tensor on which it acts. Second, when 
varying the Ricci tensor, 

5R a p = ~08g a p + V 7 V (Q %) 7 - ^V a Vp5g , (18) 

one must keep in mind that while the last two terms yield a total derivative for the Einstein- 
Hilbert term, they give a non-trivial contribution in the variation of the curvature-squared 
term. 

At the end of the day the modified vacuum Einstein's equations are given by: 

+ (Rap - \g«pR) i^r) R + { 2S (Jp) - 9^9 1& ) w 7 (^-j==^) G P s , (19) 

where the Z's are defined below. The left hand side is just as in Einstein gravity, except 
for the coefficient which is recognized as an effective, scale-dependent Newton's constant, 
(^(D) ~ Q (D£ 2 ), enforcing weaker gravity at short distances. 

The terms quadratic in curvature on the right hand side are necessary to maintain the 
Bianchi identity, V a 5S / 5g al3 = 0, guaranteed by diffeomorphism invariance of the action. 
Indeed, taking the divergence of the left hand side and commuting the covariant deriva- 
tive through Q^ 1 yields a host of terms quadratic in curvature, which precisely cancel the 
contribution from the right hand side. 

The X terms are of the form 

X a/3 = 2V a B~* 5 V h A p]5 + l -g aP V K {B^ K A jS ) + 2 (VV^) A a]s + 2 (VW^,*) B a]s 
+ 2V a A &b V^B p \ - W a B^W s A y5 - B^UA^ + A ia DB^ , (20) 

where symmetrization under (a, (3) indices is implicit everywhere. The tensors A and B to 
be substituted in (20) are in turn given by 



4(1) - 

A aj3 — 


( □ ) R ^ ' 


n af3 




4(2) _ 
— 


L 4 HQ S 1 R a/ 3 ; 


B {2) 

n af3 


- n G aP 


4(3) _ 
A aj3 — 




n a/3 


= Q s G af 3 . 



Finally, since we are interested in the cosmological evolution, we can focus on the trace 
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of (19): 

= (DL 2 ) R - 2K# (^^) Rap + (^^) G a , 

+ JW+ f ds{lW+lW) . (22) 
Jo 

One might worry that taking the trace can yield solutions which implicitly hinge on a non- 
vanishing trace-free stress tensor, such as radiation. This is of no consequence here since 
the solution of interest is self-accelerating — a radiation component would rapidly redshift 
away. 

Given the level of complexity of the above equations of motion, it is not a priori obvious 
how to extract approximate solutions. In the remainder of the section we therefore provide 
three alternative ways of deriving our self-accelerating background: i) by direct substitution 
of the cosmological ansatz in a slow-evolution approximation (Sec. II B); ii) by assuming 
that the Ricci tensor is that of an Einstein space plus a small correction (Sec. II C); Hi) by 
substituting a metric ansatz which is de Sitter plus perturbation (Sec. II D). 



B. Slow-Evolution Approximation 

Our key approximation assumes that the universe has slowly-evolving Hubble parameter: 
\H/H 2 \ <C 1, allowing us to neglect all terms of order H 2 , H, etc. This greatly simplifies 
the equations of motion. 

Starting with the leading term in (22), we have 

Q- l R= (l + D 2 L 4 + ...)i?^ 12H 2 + QH, (23) 

where the last step follows from D 2 R ~ O (jJ 2 ,H^j. Similarly, the second term is straight- 
forward to compute: 

-2R^ (^=^) Re* « {8H 2 L 2 ) 2 H. (24) 

The third term, however, is more subtle. Taylor-expanding the form factor and com- 
puting each term explicitly, we note the following structure: 2V a V /3 (DL 2 ) 3 L 2 G a p 
-(9/2) (8H 2 L 2 ) j+1 H, for all j > 1. Resumming thus gives 

2V-V 3 (^r) G ^ « ~\ (S-\SH 2 L 2 ) -1)H. (25) 
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Finally, some more algebra reveals that the jW's in (22) combine to give the following 
contribution at order H: 



J« + f ds {1^ + J< 3 >) » - — (8H 2 L 2 ) 2 H, 
Jo 2 



(26) 



Combining (23)-(26), we obtain our modified cosmological equation 



w 12/f 2 + 6ff + — (8H 2 L 2 ) 2 if + - {Q' 1 (8H 2 L 2 ) - l) if . (27) 

As a quick check, note that in the limit iff — > we recover the standard vacuum equation: 
R = 0. However we are interested in the opposite limit, i.e., HL ^> 1, in which the form 
factor Q~ x is exponentially large, and (27) reduces to 

- -?<? (8if 2 f 2 ) . (28) 

The right hand side is indeed <C 1, confirming the validity of the slow-evolution approxima- 
tion. This describes a self-inflating universe with effective slow-roll parameter 

e= - JL „ *g ( 8 if 2 f 2 ) « 1 . (29) 

As mentioned earlier, this general definition of e in terms of H reduces to the usual e 
parameter in the case of scalar-driven inflation. 



C. Perturbing Ricci 

Our solution suggests an alternative derivation in terms of a perturbative expansion about 
de Sitter space: 

R a/3 = 3H 2 g afS + r a/3 , (30) 

where ff is constant, and r is a small correction characterizing deviations from de Sitter. To 
be precise, our approximation boils down to assuming a nearly constant Hubble parameter, 
which we can write as H(t) — H + h(t), with h <C H . The components of the full Ricci 
tensor at linear order in h are then given by 

R° = 3 (if 2 + if) « 3if 2 + QHh + 3h; 

R. = (3H 2 + if) « (3if 2 + 6ff/i + /i) <f ^ , (31) 
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from which we can read off the components of r: 

r° = 6Hh + 3h + (h 2 ) 

r l . = (&Hh + fij S i j + O (h 2 ) . (32) 

By adding total derivative terms to (15), which therefore do not affect the variational 
principle, we can rewrite our action in the equivalent form 

«a f d , x ^ l R + _ r A . (33 ) 



2 J v a \ \ T ) \ □ 

Since r is assumed small, it suffices to keep terms at most linear in r when varying the 
action. After taking the trace, the equation of motion is approximately given by 

k, — 12H 2 - Q~ x r + 2V Q V /3 ^ - g ~ r af3 . (34) 

Substituting r given in (32) and, as before, neglecting terms of order h 2 , h, etc., we obtain 

h » -\h 2 G (8H 2 L 2 ) , (35) 

which is easily seen to agree with (28) by making the substitutions h = H and H — > H(t). 

D. Newton's Version 

A third method to obtain our solution, which will prove useful when calculating the 
spectrum of density perturbations, starts out with a nearly de Sitter metric ansatz, 

ds 2 = _ e 2<P(t) dt 2 + e -24>(t) e 2Ht dx 2 ^ (36) 

where is assumed small. At linear order, exp(20) 1 + 20, and we recognize as a 
time-dependent Newtonian potential. 

Substituting into the first of Eqs. (33) and truncating at quadratic order, we obtain the 
Lagrangian density 

C 



M 2 X 



= \2H 2 <P (1 + 0(0)) - 30£~ 1 D0 - 20 ^ - l - (D 2 + 18H 2 D + 36H A ) 
+ 2 (V Q V^0) (^-^) V«V^ + 0(0 3 ) , (37) 
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where we have dropped an irrelevant constant term. In other words, we treat as a scalar 
field evolving in a background pure de Sitter metric with Hubble constant H. The resulting 
equation of motion is then 

_^-i D0 _ (^j==-^ (D 2 + 18H 2 D + 36tf 4 ) 

+ V«V' 3 (^-^) V « V ^ = - 3 ^ 2 (! + W)) • (38) 

Once again we make a "slow-roll" approximation, keeping terms at most of order 0. In 
particular, this means we can neglect D 2 and higher derivative terms. Moreover, based on 
our previous derivations we anticipate that this is a consistent approximation when HL ^ 1, 
which we therefore assume from the outset. In this limit the dominant contribution from 
the left hand side comes from the last term: 

V^V 3 (^=^) V a V,0 « - (8^ 2 L 2 ) . (39) 

Thus, to leading order the equation for reduces to 

3H<j) « 8# 2 £ (8H 2 L 2 ) , (40) 

which is reminiscent of the slow-roll evolution of the inflaton: 3H<p ~ — V^,. It will indeed 
be useful to think of the Newtonian potential as an inflaton in this context, especially when 
studying density perturbations. To make contact with our previous derivations, note that the 
Hubble parameter of Sec. II B is related through via time-reparametrization: dt — > e^dt, 
and thus H —Htfi. Substituting in (40) and dropping the bars yields (28). 

III. EXITING THE SELF-INFLATIONARY PHASE AND REHEATING 

Since Q (8H 2 L 2 ) is a decreasing function of H, while H is negative, eventually the rate of 
expansion reaches \H / H 2 \ ~ 0(1), which signals the end of the self-inflationary phase. This 
occurs when HL ~ 0(\). At that point, the slow-evolution approximation breaks down, 
and one must revert to solving the full equation to describe the cosmological evolution. This 
daunting task requires numerical work which is beyond the scope of this paper. Nevertheless, 
if H remains negative throughout — a mild assumption — then within a few Hubble times 
the Hubble parameter will drop to HL « 1. In this regime, £/—>!, and (15) reduces to 
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Einstein gravity. Thus our assumption is that the universe makes a swift transition around 
HL ~ 0(1) from the self- inflationary phase to a decelerating phase described by the usual 
Friedmann equation. 

Exiting inflation is of course not the end of the story — the scenario must also account 
for reheating the universe. Since there is no inflaton to couple directly to matter fields in 
this case, the dominant reheating mechanism relies on gravitational particle production. As 
pointed out years ago by Ford [23] and others [24], the rapid change in the gravitational field 
at the end of inflation can excite matter fields which are not conformally coupled to gravity. 
(Particles of conformally invariant fields are not produced since the backround is conformally 
flat.) This results in a reheat temperature of order H — the Hawking temperature during 
inflation — , which in our case is of order 

^reheat ~ L _1 . (41) 

Of course this mechanism requires a sharp exit from inflation. A lingering accelerating phase 
would dilute the gravitationally produced particles, leaving the post-inflationary universe in 
an unacceptably cold state. This pitfall can be avoided in our case by choosing sufficiently 
steep form factors, such as our fiducial exponential example (4). 

Nucleosynthesis constrains the reheat temperature to be at least 10 MeV, corresponding 
to an upper bound on the scale of new gravitational physics: L ^ 100 fm. We will later 
see that fixing the amplitude of density perturbations imposes a much tighter constraint of 
L ~ 10~ 22 cm. 

IV. VALIDITY OF EFFECTIVE THEORY 

At first sight our action (15) looks fine-tuned from the point of view of effective field 
theory: it is at most quadratic in Ricci, derivatives appear in a special exponential form, 
and so forth. Contrary to this expectation, in this section we argue that our self-inflating 
solution is in fact robust to a host of generic corrections to (15). 

First consider adding terms to (15) involving covariant derivatives acting on curvature 
tensors, such as 

R aP W a W p R. (42) 
(This is an interesting example because of its different tensor structure than the form factor 
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in (15).) Because of the derivatives, such terms vanish in the de Sitter limit and thus can 
at best be of order H. Note that most terms will be higher order in H and thus negligible, 
such as R af3 V a RV pR ~ H 4 H 2 . The above example gives 

R aP V a V p R « -21QH 4 H . (43) 

But even so, this is still negligible compared to the exp (64if 4 L 4 ) H contribution from the 
form factor. Thus generic corrections at order H have a much smaller coefficient than that 
of the exponential form factor. 

More worrisome are derivative-free terms, such as 

R 2 L 2 ~ H A L 2 . (44) 

Since HL ^ 1 for our self-inflating solution, this is at least as important as the leading 
Einstein-Hilbert term. To see how it affects things, consider replacing the H 2 term in (27) 
with something of order H 4 L 2 , thereby modifying (28) to something of the form 

^ ~ H 2 L 2 Q {8H 2 L 2 ) . (45) 

While the H 2 L 2 factor makes the right hand side larger, nevertheless the exponential factor 
wins out and preserves the slow evolution approximation. 

To summarize, corrections to the effective action that involve derivatives acting on curva- 
ture vanish in the de Sitter limit and thus are at best of order H. Even so their coefficient is 
generically negligible compared to the exponentially large form factor. Derivative-free cor- 
rections are of the same order, or larger, than the Einstein-Hilbert term and thus cannot be 
neglected. However their influence in the equation of motion can effectively be absorbed in a 
"renormalization" of the form factor Q. Generically this renormalized form factor still van- 
ishes in the UV, preserving the slow-roll approximation. One potentially important caveat 
to this analysis is whether the standard rules of effective field theory apply in this case. 
After all our original fading gravity action displays an infinite number of time derivatives 
and is therefore genuinely non-local. 

Be that as it may, a remarkable property of our effective action is that it becomes in- 
creasingly robust with increasing HL. That is, the higher the scale of inflation compared to 
the scale of new physics, the closer our solution is to pure de Sitter, and the less sensitive it 
becomes to generic corrections. 
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V. DENSITY PERTURBATIONS 

The growth of density fluctuations is most easily studied by perturbing (36) in Newtonian 
gauge, 

rfs 2 = _ e m) (i + 2 <S>(x, t)) dt 2 + e- 24,{£) (1 - 2$ Or, t)) e 2St dx 2 . (46) 

Implicit in this ansatz is our neglecting anisotropic stress, which greatly simplifies the calcu- 
lation. Although we have not checked this in detail, this approximation is likely consistent 
since we are restricting the analysis to vacuum equations of motion only 

Postponing a careful analysis of this issue to future work, the virtue of neglecting 
anisotropic stress is that $(x, t) is now recognized as a perturbation in (pit): <p — > 
<p(t) + $(x, t). Hence the quadratic lagrangian follows immediately from (37), 

^ = S^n® -2$( 5 '"M (D 2 + \m 2 U + 36# 4 ) $ 

m 2 x \ □ ; v > 

+ 2 (V Q V /3 $) V Q V^$ + H 2 ($ 2 ) , (47) 

where the precise coefficient of the last term is not important. The corresponding equation 
of motion is 

-3^ _1 D$ - 2 (D 2 + 18H 2 D + 3QH 4 ) $ 

+2V a V /3 1 D ~ 1 ^ V Q V^$ = i? 2 ($) . (48) 

In the short-wavelength regime, fc ^> aH, we can neglect space-time curvature and set 
H to zero. Moreover, covariant derivatives commute in this limit. Thus the quadratic 
lagrangian reduces to 

£quad -> SM^g-'n® , (49) 

corresponding to Q~ l (nL 2 ) = 0, where □ now stands for the d'Alembertian in flat 
space. The canonically-normalized field variable, <j> c (x,t) = v / 6M P i^~ 1 / 2 (D)<I>, therefore 
satisfies the usual Klein-Gordon equation. Since the form factor is unity on-shell, however, 
the adiabatic vacuum corresponds to the initial amplitude 

^^br" (50) 

In the long-wavelength limit, k <C aH : the behavior is analogous to that of <p{t) described 
above since <3> just represents a time-dependent shift of the background solution. Making a 
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similar slow-roll approximation as in (39), we obtain 3H& k ~ H 2 Q (8f/" 2 L 2 ) $ fe , and thus 

$ k (t) w const. (51) 

for modes well outside the horizon. 

Note that the short and long wavelength asymptotic solutions given above match those 
of inflationary scalar perturbations: the Newtonian potential oscillates as in flat space when 
the modes are well-inside the horizon, and asymptotes to a constant when the modes exit 
the horizon. Thus, taking into account the factor of \^6Mp\ difference in normalization from 
a canonical scalar field, the spectrum on super-horizon scales is given by 

A$ fc ^— =-?-. (52) 

Alternatively we can derive (52) using the horizon-crossing approximation, which corre- 
sponds to matching (50) and the long-wavelength $ fe ^ const, solution at k = aH. In 
standard inflationary calculations this is an accurate approximation whenever H is slowly 
varying [36], which is the case here. 

A useful gauge-invariant variable to track is the curvature perturbation on comoving 
slices, denoted as usual by ( [37, 38], since it is conserved on super- horizon scales, barring 
entropy perturbations [39]. It is related to the Newtonian potential by 

H 2 ( $ H \ 
< = - T (* + ---^. (53) 

Since $ ^ const, on large scales, and H / H 2 <C 1, substitution of (52) yields 

1 H 

AC k ^^= , (54) 

where e was introduced in (29). 

As mentioned earlier, this differs by a factor of e 1 / 2 from the usual answer in scalar-driven 
inflation, tracing back to a difference in normalization of $ in the two cases. This suggests 
that the tensor-to-scalar ratio is proportional to e 2 in this case, leading to a suppressed level 
of gravity waves. We will see in Sec. VI, however, that unlike inflation scale invariance for 
tensors is not generic in fading gravity models, although certainly possible. In particular, a 
sufficiently red tilt can easily overcome the e 2 suppression and yield a measurable level of 
gravitational waves on large scales. 



19 



A. Fixing the Amplitude 



In this section we match our expression for ( with the observed large-scale amplitude of 
1CT 5 from observations of the microwave background temperature anisotropy. To do so we 
need an expression for H(N) and e(N) = —d\nH/dN, where N is defined as the number of 
e-folds before the end of inflation: 



(55) 



For the fiducial form factor (4), the results of numerically integrating (28) are shown in 
Figs. 1 and 2. A useful fitting formula for the slow-evolution parameter is 

1 



25N 



(56) 



also shown in Fig. 2 for comparison. This agrees with the exact result to within 20-40% 
except for the last few e-folds. Substituting in (54) and using the fact that H ~ L during 
inflation, the amplitude of perturbations on the largest scales must therefore satisfy 



-5 



pi 



(57) 



where N obs stands for the e-fold mark when the observable range of modes exits the horizon. 
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FIG. 1: Evolution of the Hubble parameter, H(N), during self-inflation for Q = exp (— D 2 L 4 ). 
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FIG. 2: Evolution of the effective slow-evolution parameter, e(N) = —dlnH/dN, for the same 
choice of Q as before (solid line). The fitting formula (56), eg t (-/V) = 1/25 AT", is shown for comparison 
(dashed line). 

All we need at this point is a relation between A/ obs and L. The former can be expressed 

as 

111 ( ^cnd-^cnd j _|_ J n | ^reheat -^reheat J _|_ J n | ^eq-^cq j (58) 

\ ^reheat -^reheat / \ ^eq-^eq 

where "eq" stands for matter-radiation equality, "0" stands for today, and other subscripts 
are self-explanatory. The biggest unknown is the leading term since it depends on the 
detailed evolution between the end of self-inflation and reheating. For simplicity let us 
assume that the background evolves in this intervening phase as a(t) ~ t 1 ^ 3 , as if dominated 
by the kinetic energy of a scalar field: 

In I ° cnd ^ cnd ] _ 2 ln ( ^end \ _ 2 } n | ^reheat 
\Oreheat -^reheat/ 3 \ -^reheat / 3 \ M P1 

In the last step we have used the fact that reheating marks the onset of the radiation- 
dominated epoch, and thus H rc ^ CSut ~ T r ^ heat /Mpi ~ H^ nd /Mp\. It is easy to check that 
the final result is insensitive to the details of the intervening phase — it only changes the 
coefficient in (59) by order one corrections. 

Coming back to (58), z COl 3000 implies that the last term in (58) is ps 4. More- 
over, since the universe is radiation-dominated between reheating and equality, we have 
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ln(a re heat-ffreheat/aeq-ff C q) ~ ln(T rcheat /T cq ) . Putting everything together we obtain: 

N ohs « 4 + In (^j + i In « 65 - ± In (LM P1 ) , (60) 

where in the last step we have used T rehcat ~ L~ x . 

We can now substitute this expression for N ohs into (57) and obtain a relation for L. Thus 
the observed large-scale amplitude of density perturbations fixes the scale of new physics to 

L~10 10 GeV, (61) 

which, unfortunately, lies well beyond the reach of laboratory tests of the gravitational 
inverse square law. Despite the relatively large number of e-folds required to explain the 
observed homogeneity and isotropy of our universe, the reheating temperature of 10 10 GeV is 
moderately low, owing to the inefficiency of the reheating mechanism. It should be stressed 
that (61) applies to our exponential fiducial form factor — more general form factors can 
yield higher reheating temperature. 



B. Scalar Spectral Index 

The power spectrum is given by P(k) ~ ( 2 ~ H 2 /e 2 , where evaluation at horizon exit 
is understood. Since e and H are time-varying during inflation, this leads to a spectral tilt 
defined as usual by 

_ dlnPjk) _ f dlnH _ dlne \ 
Us ~ 1= dink ~ 2 \dtol~dfok) ■ { ] 

The approximate correspondence between "time" and "scale" dependence is set as usual by 
the horizon-crossing condition, k = aH, which gives dink <-> dN. 

For our fiducial choice of Q, the equation of motion (29) written in terms of e implies 

dine , . /3 \ ,„„ N 
-4eln -e . (63) 



dink \8 , 

Moreover, since d In H/ dN = — e by definition, the spectral index can be expressed as 

n a -l = 2ej-l + 4]n^0j«8elnQe^ . (64) 

Numerically, we can substitute iV 60 in the fitting formula (56) to obtain e ~ 10~ 3 , 
corresponding to n s 0.96. This is remarkably close to the best fit value from the WMAP 3- 
year data [25] . Of course more general form factors will result in a spread of allowed n s — this 
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will be explored in subsequent work. However the indications are that the scalar spectrum 
in our model is degenerate with the generic single-field inflationary prediction [21, 22]: a 
red spectrum with a few percent deviation from scale invariance. Indeed, the scalar spectral 
tilt is unambiguously red if Q is monotonic. Since e increases with time in this case, while 
H < 0, it follows that H/e is a decreasing function of time, resulting in smaller amplitude 
for shorter wavelength modes. 

VI. PRIMORDIAL GRAVITATIONAL WAVE SPECTRUM 

The observational predictions of the model are so far qualitatively degenerate with the 
simplest models of scalar-driven inflation: nearly scale-invariant, adiabatic fluctuations, with 
a slightly red spectral tilt. A promising observable for distinguishability are tensor pertur- 
bations or primordial gravitational waves. 

We perturb the metric as 

9a0 = 9a(3 + Kp , (65) 

where for the purpose of this calculation the background metric g is taken to be exact de 
Sitter: g a pdx a dx> 3 = -dt 2 +e 2Ht dx 2 . The non-vanishing components for tensor perturbations 
are purely spatial, h° a = 0, and satisfy the usual transverse, traceless conditions: h\ = 0, 
dih'j = 0. The latter consist of four conditions on a symmetric tensor in 3 dimensions, 
leaving two degrees of freedom — the two polarization states of the graviton. 
Our starting point is the action written in the form (33) 

S=^J d'x^Tg \r + (G a/3 + 3H 2 g^) (^^) {R*fi ~ 3^) } , (66) 

where we have dropped the bars for simplicity. Substituting (65) and dropping terms that 
are independent of h, we obtain a quadratic action for tensor modes 

Censor = ^ J dV=^ (□ " 8^) {l + (^^j ( D " } V " ( 67 ) 

The constant piece in the kinetic term indicates that the usual graviton mode has a positive 
mass squared in this background. Already this is a significant difference from scalar-driven 
inflation where the graviton action describes two massless degrees of freedom. 

The above propagator is quite complicated and requires a careful analysis of its pole struc- 
ture. To do so, consider a basis of eigenfunctions h\f for the □ operator, with dimensionless 
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momentum eigenvalues — q : 



3 w * 3 

The inverse propagator, P -1 , is then given in terms of q by 



Uti. {q) = -8H 2 q 2 ti. {q) . (68) 



8H 2 v ^ 7 I V 
which, as illustrated in Fig. 3 for our fiducial form factor (4), has three zeros, corresponding 
to poles of the propagator. Although the figure focuses on a specific choice of Q, we can 
argue that the number and locations of the poles are generic. 

The zero at — q 2 = 1 is manifest from (69) and exists for any Q. The other zeros can be 
found by setting the factor in curly brackets to zero to obtain 

1 + q 2 ~ G (-8H 2 L 2 q 2 ) = . (70) 

This is clearly satisfied near — q 2 = 1 at approximately — q 2 « 1 — Q (8H 2 L 2 ). This, combined 
with the first pole, results in a quasi double pole near — q 2 = 1, which is a robust consequence 
of the approximation Q (8H 2 L 2 ) <C 1 required for the background self-inflating solution. 
Finally, there is a third zero at small — q 2 , denoted by — q 2 mal i, which is easily seen by taking 
the log of (70): 

g s 2 mall « log (G(-8H 2 L 2 q 2 m!lll )) . (71) 
For the fiducial form factor (4), this gives 

-9Lui«64^ w0 - 12 ' (72) 
where in the last step we have used (29) and (56) with N m 60. The location of this third 
pole is in good agreement with the numerical result shown in Fig. (3). While the double 
pole near — q 2 = 1 is robust under generic choice of Q, we see that — ? 2 mall depends on the 
details of the form factor. Since it is precisely this pole which ends up determining the 
tensor spectral index, as we will see below, our model does not make a generic prediction of 
scale invariance for the gravitational wave spectrum. 

The sign of the residue for the each of the poles can be read off from the slope of 
positive slope corresponds to "right sign" residue. Hence the pole at — q 2 ~ 1 — Q {8H 2 L 2 ) 
is ghost-like. The existence of ghost-like excitations in the momentum range where we trust 
the background solution is certainly worrisome, and we will devote Sec. VI B to analyze in 
detail the consequences of this instability. 
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f- 1 (-q 2 )/8H 2 




FIG. 3: Inverse graviton propagator in de Sitter background. Although difficult to distinguish 
in the plot, V~ 1 (q 2 ) actually crosses zero at — q 2 « 1 — Q (8H 2 L 2 ) and — q 2 = 1. Therefore the 
propagator has a total of 3 poles, with the one at — q 2 « 1 — Q (8H 2 L 2 ) corresponding to a ghost. 

While we have so far focused on helicity-2 modes of the metric, the modified propagator 
in (67) shows that our gravitons are massive and therefore are expected to have 5 polariza- 
tions. To be precise, the full quadratic lagrangian in the usual gauge, V a {h a p—g a ph/2) = 0, 
is given by 



--quad 



-h** (□ - 8H 2 ) (l + (^"^) (□ " 8 ^ 2 )) ^ - ^hUh , (73) 



where h a p = h a p — (japh/A is the traceless mode. Thus only the traceless part of the 
metric perturbation acquires a modified propagator, while the action for the trace part 
is unchanged compared to general relativity. In general we therefore expect 5 propagating 
degrees of freedom for each of the 3 poles, for a total of 16 (including the trace). Interpreting 
— q 2 as an effective mass from (68), each set of 5 polarizations should form a massive spin-2 
representation. The analysis of Sec. II D suggests that one of the corresponding helicity-0 
states plays the role of an effective irmaton in our model, setting the clock for the cosmological 
evolution. A rigorous analysis of the propagating degrees of freedom in this model is currently 
in progress [40]. 
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A. Tensor Spectrum 

Putting aside the ghost issue for the time being, let us charge ahead and compute the 
primordial gravitational wave spectrum. The calculation is closely related to its scalar 
counterpart. Well-within the horizon, the wave equation reduces to its flat space form, with 
usual Bunch-Davies solution 

fefc-> — 1 r^ ^- ( 74 ) 



M Pl aV2k 

In the long wavelength limit, we have in principle three graviton excitations to keep track 
of, corresponding to the above three poles. Fortunately, however, only the pole at small — q 2 
is relevant — the two poles near — q 2 = 1 correspond to massive gravitons and therefore are 
not excited by the background expansion. The wave equation for this mode is approximately 
given by 

-UW 3 = h) + 3Hh) - 2H% » 8H 2 q 2 miill h) , (75) 
which has the growing mode solution 

h k = C k jU 1+ ^ +4 ^^ . (76) 

The behavior at intermediate wavenumber requires solving the perturbation equation 
exactly, which is beyond the scope of this paper. For our purposes it suffices to use once 
again the horizon-crossing approximation, as we did for the scalar spectrum: we treat (74) 
and (76) as exact solutions for k > aH and k < aH, respectively, and determine C k by 
matching them at k — aH. The final answer for the large-scale gravitational wave spectrum 
is then 

k 3/2 H ( -1(^1+1 (i+H 2 ma ii)-i) 

7T 

corresponding to a tensor spectral tilt 



w ~AfeU> ■ (77) 




1 + 5 (1 + 4?L.n) - 1 J P8) 

Although apparently not generic, a nearly scale invariant tensor spectrum is possible if 
— Ismail ~ 1/4. It is straightforward to show that this can be achieved with the following 
generalization of our fiducial form factor, 

g (DL 2 ) = exp {aUL 2 - D 2 L 4 ) , (79) 
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with a ~ 0(1) suitably chosen. Even in this case, however, the amplitude of gravity waves is 
well below the sensitivity levels for near-future B-mode polarization experiments. Indeed the 
scalar amplitude (54) implies a tensor-to-scalar ratio of r ~ e 2 , which is at best 10~ 4 — 1CT 3 . 

The key implication of (78), however, is that the tensor spectrum is generically not scale 
invariant — its spectral tilt is sensitive to the choice of form factor. This is a distinguishing 
feature from scalar-driven inflation, where a nearly scale-invariant tensor spectrum follows 
directly from the slow evolution of the Hubble parameter during inflation. In particular, a 
significant gravity wave amplitude on large scales requires a sufficiently red tilt, which is 
not hard to obtain. Our fiducial form factor, for instance, corresponds to — 9smaii ~ 0.1, and 
thus ht ~ —1/3, which is in fact too red to be consistent with observations. A more general 
form factor, such as (79), can do the trick. 

The bottom line is that if gravitational waves are to be observed in our model, their 
spectrum must have a relatively large red tilt, which would distinguish it from standard 
inflation. Conversely, observing gravity waves with a nearly scale invariant spectrum and 
with a tensor-to-scalar ratio at the percent level would rule out fading gravity inflationary 
models. Of course failing to detect tensor perturbations would constrain but not rule out 
either standard or fading gravity inflation. 

If futuristic CMB polarization or direct gravity wave experiments ever get down to sen- 
sitivity levels that would probe tensor-to-scalar ratios in the 1CT 4 — 1CT 3 range, then our 
model would be strongly favored by the observation of a blue tilt for the tensor spectrum. 
This outcome, which is certainly allowed by (78), would be an indisputable distinguishing 
feature from scalar driven inflation, where the spectrum is unequivocally red. 

B. Taming the Ghost 

As usual it is possible to regularize the ghost instability by invoking Lorentz-violation at 
some high cutoff scale. We explore this possibility in this subsection. In the limit HL <: 1, 
the ghost can be thought of as arising from an approximate double-pole at — q 2 = 1, or 
□ = 8H 2 . Indeed, in this regime the quadratic action (67) near □ = 8H 2 is 

Aensor ~ G' 1 {8H 2 L 2 ) (□ - Stf 2 )' h tJ . (80) 
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Incidentally, because of the Q~ x (8H 2 L 2 ) prefactor the ghost decouples in the limit HL — > oo. 
Unfortunately, as we will see in detail below, this decoupling limit leads to unacceptably 
large density fluctuations. The COBE constraint on the amplitude of 5p/p forces us to take 
the ghost by the horns. 

A ghost indicates an explosive instability of the theory. Because ghost particles carry 
negative energy, the vacuum is unstable to decay, e.g., through the process shown in Fig. 4 
in which a ghost graviton and two matter particles are spontaneously created from the 
vacuum. It is well known that the rate for this decay is formally infinite due to divergent 
phase space factors, even with a Lorentz invariant cut-off. The reason is simple: the energy 
carried by the ghost can be made arbitrarily large while keeping p 2 small, therefore satisfying 
some Lorentz invariant bound. 




FIG. 4: Vacuum decay to ghost graviton and matter particles. 

The decay rate can be regularized by invoking explicit Lorentz breaking at some scale 

Ali (measured in the cosmic rest frame) which sets an upper bound on the magnitude of the 

ghost energy [32]. Ignoring irrelevant prefactors, the decay rate per unit volume is given by 

A 6 

r ^ s vv> ~ JffG (8H 2 L 2 ) , (81) 
where the decoupling factor has trickled down from the propagator (80). 

Our self-inflating solution remains valid if the ghost energy density produced in a Hubble 
time, 

(82) 

is less than the background effective energy density given by 3H 2 M^ X : 
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This condition ensures that the rate of ghost production is sufficiently slow so as not to 
drown the inflating background. Self-inflation keeps going on in this case, continuously 
diluting the ghost and matter particles created spontaneously from the vacuum. 

It seems a priori that this bound can be trivially satisfied even for Ali ~ Mpi by imposing 
a sufficiently high scale for inflation. However such a limit also makes the amplitude of 
perturbations too large. To see this explicitly let us substitute the expression for Sp/p 
using (54) and (29): 



Thus, unfortunately this has the decoupling factor with an inverse power. Substituting 
5p/p ~ 1CT 5 and the value of Q (8H 2 L 2 ) for the observable range of modes — see Sec. VA 
— yields a lower bound on the scale of Lorentz symmetry breaking: 



Since COBE also constrains L ~ 10 10 GeV — see (61) — , this leaves us with a wide range 
of energy scales within which our effective description is valid. 

Once inflation ends the ghost further decouples from all other fields. By construction at 
the perturbative level our theory has two propagating degrees of freedom about flat space, 
which are just the usual polarizations of the graviton. Thus in the limit HL — > the 
ghost becomes non-perturbative and therefore decouples in weak field. This situation is 
in stark contrast with ghosts in other well-known higher curvature theories, such as R + 
CiR 2 + c 2 R a i3R a P [33]. In that case one generically finds an extra massive spin-2 field, which 
remains in the perturbative spectrum for arbitrarily large momenta. 

VII. GHOST STORIES AND THE OSTROGRADSKI INSTABILITY 

The ghost instability described in the previous section is a specific manifestation of a gen- 
eral pathology afflicting higher derivative theories known as the Ostrogradski instability [18]. 
The basic statement is that the Hamiltonian of such theories generically depends linearly on 
all but one conjugate momentum variables and, hence, is unbounded from below. See [19] 
for a nice exposition. The result applies to theories with finite number of time derivatives, 
but will generically survive the infinite-derivative limit for theories that are analytic in the 




(84) 



A L t < 10 15 GeV . 



(85) 
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derivative operator. An example of the latter which will prove useful is 



e DL2 D + m 




(86) 



Such analytic form factors are ubiquitous in string theory. They appear in the truncated 
tachyon action in bosonic string field theory and in p-adic string theory. More generally such 
a structure is expected from the finite UV behavior. And indeed, it remains mysterious how 
string theory avoids the aforementioned instabilities [19]. The footprints of Ostrogradskian 
instabilities are found in time-dependent solutions of p-adic actions, where the scalar zooms 
by the tachyon vacuum and undergoes oscillations of ever-growing amplitude [41]. In bosonic 
open string field theory, however, it has been argued recently that these are artifacts of the 
level truncation approximation, and that the full string field theory, albeit non-local in space, 
is local in lightcone time [42] . See [43] for an analysis of degrees of freedom in p-adic and 
string field theory. 

A related worry for higher-derivative theories is the lack of a well-defined Cauchy problem. 
To get a unique solution from an equation of motion with infinite number of time derivatives 
naively requires an infinite amount of initial data, which therefore amounts to arbitrarily 
specifying the solution over any finite time interval. 

These results afflict a wide class of non-local theories, including our model of fading 
gravity. However we would like to argue that the pathologies are of a mild form, at least in the 
regimes of interest for the cosmological analysis. Such considerations lead us to believe that 
in the neighborhood of the self-inflating solution there should exist an equivalent description 
of the theory which is local time and thus has a well-defined initial value formulation. 



An important exception to Ostrogradski's results is a non-local theory which can be made 
local by suitable invertible (non-local) field redefinitions. For instance, 



A. Local Perturbative Locality 




(87) 



is manifestly local in terms of a new field variable <p = e DZ/2 / 2 0. However, adding local mass 



term and interactions, 



OO 




n=2 
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makes the non-locality real and introduces Ostrogradskian pathologies. Nevertheless, as long 
as we restrict ourselves to perturbative analysis, with £ as the unperturbed Lagrangian, the 
theory has a well-defined initial value formulation: all perturbative solutions only depend 
on the unperturbed solution and its time-derivative evaluated at some initial time. This can 
be made manifest by constructing an alternative theory which is local in time (albeit still 
non-local in space) and reproduces all perturbative solutions of the original theory. Such a 
construction is called localization. 

The point, of course, is that the runaway solutions of the full theory are genuinely non- 
perturbative. Perturbation theory effectively projects out these dangerous solutions and 
focuses on a subspace of solutions which is stable and spanned by two pieces of initial data. 
Theories of the form (88) are therefore said to be perturbatively local. While the perturbative 
regime circumvents the pathologies of the full theory, it is not completely worry-free since its 
solutions display acausality [19], in the sense that corrections to the unperturbed solution 
at a given space-time point have support outside the past light cone. 

It is clear that the weak-field limit of our model about flat space satisfies perturbative 
locality: 

Afield = ^ j d'x | (h* - Q- 1 (DL 2 ) Uhap + 0(h 3 ) } . (89) 

In particular, one can make an invertible metric redefinition to render the kinetic term canon- 
ical, thereby pushing the non-local factors in the interaction terms. Even in more general 
backgrounds, the weak-field, flat space approximation is always valid within a sufficiently 
small region around any point. Thus our theory can be said to satisfy local perturbative local- 
ity — the perturbative UV limit is free of the Ostrogradski instability and has a well-defined 
initial value formulation. 

Beyond weak-field, however, the perturbed action about general backgrounds is not gener- 
ically perturbatively local. Consider for instance our quadratic action (67) for transverse, 
traceless perturbations about the self-inflating solution. Because of the extra poles in the 
propagator, in this case the non-locality cannot be eliminated through an invertible field 
redefinition. This is akin to (86). The Ostrogradski instability is signaled by the existence 
of the ghost-like mode discussed in Sec. VI. 
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B. Non-Perturbative de Sitter Solution 

To summarize, our toy model action (15) for weaker gravity displays local perturbative 
locality but suffers from the non-perturbative Ostrogradskian instability Moreover the 
action is non-local in time and is therefore bereft of a Cauchy problem. The non-perturbative 
nature of the self-inflating solution makes one wonder whether it secretly exploits unwanted 
instabilities of the theory. A related worry pertaining to the temporal non-locality of the 
full theory is whether the de Sitter solution requires fine-tuning an infinite number of time 
derivatives, and whether we could have instead obtained any other background evolution. 

While these questions are certainly warranted, it is encouraging to note that the self- 
inflating solution appears to rely in the most minimal way possible on the higher-derivative 
nature of the theory. Indeed, despite having an infinite number of time derivatives, the 
modified cosmological equation (28) is nevertheless a second-order differential equation in 
the scale factor, just as in Einstein gravity. This suggests that it should be possible to 
construct an alternative theory which would be local in time (although likely non-local 
in space) and would reproduce all dynamics within a solution-space neighborhood of the 
self-accelerating solution. 

These considerations naturally lead us to ask whether our self-inflating scenario is de- 
generate with some (local) scalar field model. While we cannot completely rule out this 
possibility, the allowed strong departure from scale invariance in the tensor spectrum sug- 
gests that such a scalar-field analogue is impossible. Our gravity wave spectrum can be blue. 
In Einstein gravity, on the other hand, a blue tilt for tensors requires either a contracting 
universe (as in ekpyrotic/cyclic models) or an expanding universe driven by a fluid which 
violates the null energy condition. 

The former is unlikely since contracting universes generically lead to a big crunch sin- 
gularity; although we haven't ruled out a singularity at the end of the self-inflating phase, 
everything so far points towards a graceful exit and smooth matching onto a radiation- 
dominated era. Degeneracy with null-energy violating models is also unlikely. First of all, 
the only consistent candidate known to this author is the recent proposal of [30] based on 
ghost condensation [31]. In this case, consistency of the model requires H <C H 2 , which 
yields only a mild blue tilt. Even if one could extend the validity of the model such that a 
strong blue tilt obtains, this would likely translate into large deviations from scale invari- 
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ance in the scalar spectrum as well. Degeneracy with local scalar field models is therefore 
unlikely. If a temporally-local alternative description of the self-inflating solution is possible, 
it will likely retain the spatial non-locality. 

VIII. CONCLUSIONS AND FUTURE DIRECTIONS 

In this paper we have studied the cosmology of a toy model where the gravitational 
force becomes exponentially weak at short distances. This is accomplished by modifying the 
graviton propagator by an analytic form factor which shuts off in the UV. In this context we 
found novel inflationary solutions that do not rely explicitly on scalar fields or other form of 
stress energy. Instead, inflation results from the modified vacuum equations of motion. Of 
course there are new degrees of freedom associated with the modified propagator, in partic- 
ular a scalar mode which acts as a clock for the inflationary evolution. What is intriguing, 
however, is that this effective inflaton effectively decouples and therefore disappears from 
the perturbative spectrum in the flat space limit. 

A novel inflationary mechanism is interesting because it opens the door to a new realm of 
"generic" predictions. Moreover, what may appear as fine-tuning in scalar-driven inflation 
could be more natural in the fading gravity context, and vice versa. 

A remarkable feature is that de Sitter space is an increasingly accurate solution to the 
fading gravity equations as HL — * oo. Generic corrections to the effective action become 
less and less relevant for the background solution in this limit. Usually the exact opposite 
happens: corrections to the effective action become more important above the scale of 
new physics. The stability of our solution originates from the special property that the 
Ricci curvature tensor for pure de Sitter vanishes when acted upon by covariant derivative 
operators. Thus, pure de Sitter is oblivious to the form factor, but small deviations from it 
are highly sensitive to its action in the limit of large HL. 

This work was mostly exploratory in nature, and many important issues must be resolved 
for the scenario to be viable: 

• The details of the transition from self-inflation to a radiation-dominated universe must 
be worked out in detail. There are compelling reasons to believe such a transition is 
possible: self-inflationary dynamics occur with \H\/H 2 increasing, so that the slow- 
evolution approximation eventually breaks down; radiation is produced through par- 
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tide production; as HL — > the theory reduces to Einstein gravity. While these are 
compelling hints, it is nevertheless conceivable that the self-inflationary phase and the 
normal radiation-dominated evolution actually lie on different branches of solutions or 
are separated by a curvature singularity. This is currently under investigation [40]. 

• The perturbation analysis around the self-inflating solution reveals a massive tensor 
mode which is ghost-like, indicating a violent instability. We have used standard argu- 
ments invoking Lorentz symmetry breaking at high energy to regularize the instability 
and render our solution viable. A better understanding of the origin of this ghost 
might suggest an improved, ghost-free version of the model. 

• A distinguishing prediction is that the spectrum of gravity waves is generically not 
scale invariant. In particular a relatively strong red tilt is possible, without jeopar- 
dizing the scale invariance of the scalar spectrum. It would be interesting to study in 
greater detail the observational constraints, as well as the implications for microwave 
polarization experiments and gravitational wave detectors. 

• From a phenomenological perspective, other potential observable signatures should be 
investigated, such as the predicted level of non-gaussianity. 
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